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Abstract 

In this paper we present a versatile method for the investigation of interaction net- 
works and show how to use it to assess effects of indirect interactions and feedback 
loops. The method allows to evaluate the impact of optimization measures or fail- 
ures on the system. Here, we will apply it to the investigation of catastrophes, in 
particular to the temporal development of disasters (catastrophe dynamics). The 
mathematical methods are related to the master equation, which allows the appli- 
cation of well-known solution methods. We will also indicate connections of disaster 
management with excitable media and supply networks. This facilitates to study 
the effects of measures taken by the emergency management or the local operation 
units. With a fictions, but more or less realistic example of a spreading epidemic 
disease or a wave of influenza, we illustrate how this method can, in principle, pro- 
vide decision support to the emergency management during such a disaster. Similar 
considerations may help to assess measures to fight the SARS epidemics, although 
immunization is presently not possible. 

Key words: Master equation, interaction network, excitable media, supply chain 
management, robustness of graphs, causality network, catastrophe dynamics, 
disaster preparedness 



1 Introduction 



Natural disasters [1] liave occured since earliest times, and despite the devel- 
opment of science and technology, they still cause many victims each year. 
One reason for this is the increased world population. Nowadays, more people 
can be affected by a disaster than in former centuries. Apart from this, more 
people affect their environment, which challenges its stability and triggers 
disasters such as famines. 
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One common feature of many disastrous events is the so-called domino or 
avalanche effect. It means, that one critical situation triggers another one and 
so one, so that the situation worsens even more. One famous example is a 
mountain slide that fell into a lake and caused very high waves (Vajont, Italy 
1963 [2]). Other examples are fires and failures of water and electricity supply 
caused by earthquakes or, on a larger timescale, a disease harming the economy 
and social life of the affected area, which leaves the people and country even 
poorer and the medical system less effective, causing further fatalities (e.g. 
AIDS in Central Afrika [3] or the plague in former ages). 

Therefore, a great effort is necessary to prevent the emergence of disasters 
(which is not always possible) and to improve the management of catastrophes. 
Physics and other natural scientists have traditionally contributed a lot to 
understanding the laws behind catastrophes. For example, we mention the 
extensive work on forest fires [4] and earthquakes [5] , which relate to concepts 
of self-organized criticality used to describe avalanche effects [6]. But there 
is also work on fioodings [7], landslides [8], and volcanos [9]. A considerably 
attention has also been devoted to epidemics [10]. 

Some of the physical disciplines involved into the study of these subjects are 
the theory of catastrophes and bifurcations [11], non-equilibrium phase transi- 
tions [12], self-organized criticality and scaling laws [6], percolation theory [13], 
the statistical physics of networks [14] and extreme events [15], stochastic pro- 
cesses [16] and noise-induced transitions [17], but mechanics, fluid-dynamics, 
and other fields play an important role as well [18]. 

In this paper, we like to develop a flexible semi-quantitative method allowing 

• to assess the suitability of alternative measures of emergency management, 
i.e. to give decision support, 

• to estimate the temporary development of catastrophes, and 

• to give hints when to take certain actions in an anticipatory way. 

For that purpose, it is necessary to take into account all factors which are 
relevant during the catastrophe and all direct and indirect interactions between 
them. This method is in the tradition of systems theory [19]. It extends the 
concept of the causality diagram in Sec. 2, while a dynamical generalization 
is developed in Sec. 4. 

In the next section, we start with a static analysis of interaction networks. Sec- 
tion 3 is intended to illustrate its usefulness for disaster management. There, 
we will calculate in a semi-fictious example the effects of measures taken to 
combat an epidemic catastrophe. Section 4 contains the extension to a dy- 
namic description, which is connected to the discrete master equation [16]. It 
allows to determine the probability and order of events as well as their most 
likely occurence in time. In Sees. 5 and 7, we will develop a dynamical model 
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of disaster management, which specifies some parameters in the master equa- 
tion. It relates to models of excitable media [20] and supply networks [21,22]. 
Some analytical results for this model are presented in Sec. 6, while Sec. 8 
summarizes our results and closes with an outlook. 



2 Assessment of Interaction Networks 

In this section, we want to develop a simple method to refiect the approximate 
influence of different factors or sectors on each other. Such factors may, for 
example, be energy supply, public transport, or medical support. In principle, 
it is a long list of variables i, which may play a role for the problem under 
consideration. If we represent the influence of factor j on factor i by A^j, we 
may summarize these (directed) influences by a matrix A = (Aij). However, 
in practical applications, one faces the following problems: 

(1) The number of possible interactions grows quadratically with the num- 
ber of variables or factors i. It is, therefore, difficult to measure or even 
estimate all the influences Aij. 

(2) While it appears feasible to determine the direct influence Mij of one 
variable j on another one i, it is hard or almost impossible to estimate 
the indirect influences over various nodes of the graph, which enter into 
Aij as well. However, feedback loops may have an important effect and 
may neutralize or even overcompensate the direct inffuences. 

Problem (i) can be partially resolved by clustering similar variables and se- 
lecting a representative one for each cluster. The remaining set of variables 
should contain the main explanatory variables. Systematic statistical methods 
for such a procedure are, in principle, available, but intuition may be a good 
guide, when the quantitative data required for the clustering of variables are 
missing. 

Problem (ii) may be addressed by estimating the indirect inffuences due to 
feedback loops based on the direct inffucences M^j, which can be summarized 
by a matrix M = {Mij). One may use a formula such as 



but as this converges only for small enough values of r, we will instead use 
the formula 



1 CXD 1 CX) OO 




(1) 



A = A^ = - 



1 ^ t^M!" 1 



[exp(TM) - 1] , 



(2) 



k=l 
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where 1 denotes the unity matrix. The expression reflects all influences 
over k — 1 nodes and k links, i.e. k = 1 corresponds to direct influences, 
A; = 2 to feedback loops with one intermediate node, = 3 to feedback loops 
with two intermediate nodes, etc. The prefactor t'' is not only required for 
convergence, but with r < 1, it also allows to reflect that indirect interactions 
often become weaker, the more edges (nodes) are in between. 

A further simplification can be reached by restricting to a few discrete values 
to characterize the influences. We may, for example, restrict ourselves to 



where M^j = ±3 means an extreme positive or negative influence, Mij = ±2 
represents a strong influence, M^j = ±1 a weak influence, and Mij = a 
negligible influence. Of course, a flner differentiation is possible, whenever 
necessary. (For an investigation of stylized relationships, it can also make 
sense to choose G { — 1, 0, 1}, where My = ±1 represents a strong positive 
or negative influence, then.) 

The matrix A = (^y) will be called the assessment matrix and summarizes 
all direct influences (M) and feedback effects (A — M) among the investigated 
factors. It allows conclusions about 

• the resulting strength of desireable and undesireable interactions, when feed- 
back effects are included, 

• the effect of failures of a speciflc sector (node), 

• the suitability of possible measures to reach speciflc goals or improvements, 

• the side effects of these measures on other factors. 

This will be illustrated in more detail by the example in Sec. 3. 

One open problem is the choice of the parameter r. It controls how strong the 
indirect effects contribute in comparison with the direct effects. A small value 
of T corresponds to neglecting indirect effects, i.e. 



while increasing values of r reflect a growing influence of indirect effects. This 
is often the case for catastrophes, as these are frequently related to bifurcations 
or phase transitions, to avalanches or percolation effects [5,4]. By variation of 
T, one can study different scenarios. 

Note that r may be interpreted as time coordinate: Deflning 



M,,- £{-3,-2,-1, 0,1, 2, 3}, 



(3) 



lim A^ = M , 



(4) 




(5) 
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for an arbitrary vector X, we find -'^(O) — X, 



T T 



and 



dX Xr -XO 



Prom this point of view, 

X{t) = (rA, + l)X(O) (6) 

— * — * 

describes the state of the system at time r, and Mij the changing rates. X = 
is a stationary solution and corresponds to the normal (everyday) state. An 
initial state A(0) 7^ may be interpreted as perturbation of the system by 
some (catastrophic) event. We should, however, note that the hnear system 
of equations (6) is certainly a rough description of the system dynamics. It 
is expected to hold only for small perturbations of the system state and does 
not consider damping effects due to disaster management. Such aspects will 
be considered later on (see Sees. 4 and 7), after discussion of an example 
illustrating how to apply interaction matrices to cope with catastrophes. 



3 Optimization of Interaction Networks: A Simple Example 



One advantage of our semi-quantitative approach to catastrophes is that it 
allows to estimate the impact of certain actions on the whole set of factors. 
Usually, during a disaster the responsibilities have only dissatisfactory infor- 
mation and short time to decide, so in many cases they will take into account 
only direct impacts on other factors. In the worst case, this may lead to the 
opposite than the desired result, if the feedback effects exceed the direct in- 
fluence. Therefore, it would be better to know the implications on the whole 
system. As we have argued before, all direct and indirect effects are sum- 
marized by the matrix A, which is determined from the matrix M of direct 
interactions. Different measures taken by the responsibilities are reffected by 
different matrices M. 

As an example we consider the spreading of a disease. For illustrative reasons, 
we will restrict to the discussion of five factors only: 

(1) the number of infected persons, 

(2) the quality of medical care, 

(3) the public transport, 

(4) the economic situation and 

(5) the disposal of waste. 
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These factors are not independent from each other, as illustrated by Fig. 1. 




(7) 



Figure 1. Interaction network for the example of a spreading disease discussed in 
the text. 

The corresponding matrix of the assumed direct influences among the different 
factors is 

-2+2 -1^ 
-2 +1 +2 +1 
M= -10 0+2 
-10+2 +1 
y-1 +1 +2 y 

The right choice of the sign of the direct influence M^j of factor j on factor i is 
plausible: We assume a positive sign if the factor i increases with an increase 
of factor j, while we assume a negative sign when factor i decreases with the 
growth of factor j. However, the determination of the absolute value of Mjj 
requires empirical data, expert knowledge, or experience. We have argued as 
follows: 



• A growing number of infected persons affects all other factors in a negative 
way (see first column), as these will not continue to work. That is, there 
will be problems to maintain a good economic situation, public transport, or 
the disposal of waste. Health care is affected twice, as not only the medical 
personnel may be infected, but also a higher number of patients needs to 
be treated, and capacities are limited. Therefore, we have chosen the value 
—2 in this case, but —1 for the other factors. 

• A well operating health system (second column) can reduce the number 
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of infected persons efficiently, so that we have chosen a value of —2. The 
influence of the health system on the economic situation and other factors 
was assumed to be of indirect nature, by reducing the number of ill persons. 

• Public transport (third column) contributes to a fast spreading of the in- 
fection assumed here. Therefore, we have selected a value of 2. Transport 
is also an important factor for economic prosperity (therefore the value of 
2), and it is required to get medical personnel and workers in the disposal 
sector to work (which is reflected by a value of 1). 

• The economic situation (fourth column) has a significant effect on the qual- 
ity of the health system, public transport, and disposal, so that we have 
chosen a value of 2 in each case. 

• Waste may contribute to the spreading of the disease, if it is not properly 
removed. Therefore, a good disposal system (fifth column) may reduce the 
number of infections (therefore the value of —1). It is also required for a 
functioning health system and economic production. That is, why we have 
assumed a value of 1. 

Depending on the respective situation, the concrete values of the direct in- 
fluences Mij may be somewhat different. For their speciflcation, it can also 
be helpful to check the resulting values of Aij of the overall direct plus indi- 
rect influence for their plausibility, and to compare the size of second-order 
or third-order interactions. For example, we see that the third-order feedback 
loop "number of infected persons^economic situation— s>quality of the health 
system— ^number of infected persons" is proportional to (—1) • (-1-2) • (—2) = 4. 
The same indirect influence is found for the feedback loop "number of infected 
persons— ^economic situation— > public transport— > number of infected persons" . 
Moreover, according to our assumptions, the second-order autocatalytic in- 
crease of the number of infected persons via its impact on the health system 
is four times as large as the one via its impact on the waste disposal. One 
surprising observation is that the number of infected persons is reduced via 
its impact on pubhc transport. In fact, once less busses are operated (because 
the bus drivers are ill), the spreading rate of the disease is reduced. This may 
inspire responsibilities to reduce public transport or even stop it. Later on, we 
will discuss the effect of this possible measure. 

Before, let us have a look at the resulting overall interaction matrix 



0.9 -2.2 1.3 -0.8 -1.6 



-3.4 



1.1 1.5 



3.5 2.3 




(8) 



-2.0 



0.6 2.1 



1.5 1.6 



-2.0 



0.6 1.5 



2.9 0.9 
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For its calculation, we have chosen the value r = 0.4, which will also be used 
later on to assess alternative actions to fight the spreading of the disease. 
In order to discuss a certain scenario, we will assume that Xj reflects the 
perturbation of factor j. Because of Eq. (6), the quantities 



will be used to characterize the potential response of the system in the specific 
scenario described by the perturbations Xj (and without the damping effects 
by disaster management discussed in later sections of this contribution). Here, 
5ij denotes the Kronecker function, which is 1 for i = j and otherwise. We 
will assume Xi — 1.0, as the number of infected persons is higher than normal, 
and X2 = X^ = X4 = X^ — —0.1, as the other factors are reduced by the 
spreading of the disease: 



Moreover, if we attribute a weight of Zi = 0.5 to the number of infected 
persons, a weight Z4 = 0.3 to the economic situation, and weights of Z2 ~ 
Z3 = 0.1 to the quality of the medical care and public transport, while we do 
not care about waste in our evaluation (i.e. Z^ — 0), the resulting value of 



will be used to assess the overall situation of the system. In the stationary 
(normal) system state, F would be zero. Therefore, we want to find a strategy 
which brings F close to zero. For our basic scenario, we find 



(yi,y2,n,n,n) = (l-5,-l-8,-1.0,-l-l,-l-l) and F = 1.4. (12) 



These reference values will be compared with the values for alternative sce- 
narios which correspond to different actions taken to fight the catastrophe. 

For example, let us assume to have limited stocks of vaccine for immunization. 
Should we use these to immunize 1) the transport workers, 2) the medical staff, 
or 3) the disposal workers? In the first case, we have the modified matrix 



Yi^Y.i^Aij + Sij)Xj 



(9) 



j 



{Xi, X2, X3, X4, X5) = (1.0, -0.1, -0.1, -0.1, -0.1) . 



(10) 




(11) 



0-2+20- 



-2 +1 +2 +1 



M 



0+2 



(13) 



-10+2 +1 



V 



-1 0+1+2 
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which imphes 



A 



^ 1.2 -2.3 1.4 -0.8 -1.7^ 
-3.2 1.1 1.5 3.5 2.2 
-0.5 0.1 0.8 2.4 0.5 
-1.6 0.5 2.2 1.5 1.5 

y-1.7 0.6 1.6 2.9 0.9 



(ri,F2,n,>4,n) = (1.6,-1.7,-0.5,-1.0,-1.0), and F=1.4. 
In the second case, when we immunize the medical staff, we find 



which imphes 



M 



-2 +2 -1^ 
^0+1+2 +1 
-10 0+2 

-10+2 +1 
-1 0+1+2 



0.5 -2.1 1.3 -0.7 -1.5 
-2.3 0.7 1.8 3.4 2.0 
-1.7 0.6 0.5 2.5 0.8 
-1.9 0.6 2.1 1.5 1.6 
-1.9 0.6 1.6 2.9 0.9 



(ri,r2,n,n,n) = (1.3,-1.3,-0.9,-1.1,-1.1), and F = 1.2 . 
In the third case, when the disposal workers are immunized, we expect 



M 



-2 +2 -1^ 
-2 +1 +2 +1 
-10 0+2 
-10+2 +1 
y +1 +2 y 
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which imphes 





0.6 


-2.1 1.3 


-0.7 


-1.6' 




-3.1 


1.1 1.6 


3.5 


2.2 




-1.6 


0.6 0.5 


2.5 


0.8 




-1.7 


0.6 2.1 


1.5 


1.5 




^-0.8 


0.2 1.9 


2.8 


0.6 ^ 



(20) 



(yi,y2,l3,n,n) = (l-4,-l-7,-0.9,-1.0,-0.7), and F = 1.3. (21) 



While the immunization of the pubhc transport staff has almost no effect on 
the overall situation in the system, the last two measures can improve it. We 
see that it is more effective to immunize the medical staff than the disposal 
workers, and the best would be to immunize both groups. This corresponds 
to 



and we obtain 



M 



0-2+2 -1 
^0+1+2 +1 
-10 0+2 
-10+2 +1 
0+1+2 



0.2 - 


-2.0 1.2 -0.7 


-1.5 


-1.9 


0.6 1.9 3.4 


1.9 


-1.6 


0.5 0.5 2.5 


0.8 


-1.7 


0.6 2.1 1.5 


1.5 


-0.8 


0.2 1.9 2.8 


0.6 


-1.2, 


-0.9,-1.0,- 


0.6), 



(22) 



(23) 



Other measures do not change the interactions in the system, but correspond 
to a change of the effective impact X of the catastrophe. For example, we may 
consider to reduce pubhc transport. With (8) and 

(Xi,X2,X3,X4,X5) = (1.0,-0.1,^,-0.1,-0.1), (25) 

we find 

(Fi, 1^2, 1^3, n) = (1.0,-2.4,-2.0,-1.9,-1.7) and F = 1.6. (26) 
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We see that the number of infections could, in fact, be reduced. However, the 
overall situation of the system has deteriorated, as the economic situation and 
all the other sectors were negatively affected, because many people could not 
reach their workplace. Therefore, let us consider the option to increase the 
number of disposal workers. With (8) and 

(Xi,X2,X3,X4,X5) = (1.0,-0.1,-0.1,-0.1,M), (27) 

we find 

(Fi, 1^2, 1^3, n, n) = (1.1, -1.3, -0.8, -0.8, -0.3) and F = 1.0. (28) 
In conclusion, increasing the hygienic standards can be surprisingly efficient. 

Finally, let us assume an improved waste disposal together with the immu- 
nization of both, the medical staff and the disposal workers. In that case the 
interactions of the relevant factors are characterized by matrix (22), whereas 
the starting vector is again (27). The resulting response is 

{Yi, Y2, ¥3, 14, F5) = (0.8, -0.7, -0.7, -0.6, 0.1) and F = 0.75 . (29) 

Only this combination of measures manages to actually reduce the infections 
compared to the initial state, i.e. Yi < Xi. However, we can also see that a neg- 
ative impact on the economic situation and other factors cannot be avoided. 
In any case, we can assess which measures are reasonable to take, which im- 
pact they will have on the system, and which measures need to be combined 
in order to control the spreading of the disease (or other problems in different 
scenarios) . 

The simple example in this section was chosen to illustrate the procedure 

how to assess interaction networks and potential optimization measures. In 
an on-going project, we do now investigate the interaction network among a 
large number of factors for the floodings in Germany during August 2002 and 
other catastrophes. This involves considerably more detailed and much larger 
matrices M, where nobody would be able to assess the feedback loops without 
a method such as the one proposed above. It also involves other aspects such 
as human forces fighting the catastrophe and the availability of technical or 
other equipment etc, which will be modeled in the following sections. 



4 Impact of the Interaction Network on Catastrophe Dynamics 

Before, we have used the interaction network predominantly for the static as- 
sessment of the influence of different factors on each other. We will now try 
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to extend this method step by step in a way that allows a semi- quantitative 
analysis of the time-dependence of catastrophes for the purpose of anticipa- 
tion, which helps to prepare for the next step in catastrophe management or 
prevention. We are particularly interested in the domino or avalanche effects 
of particular events such as the failure of a certain factor or sector in the in- 
teraction network. We will assume that this faihire spreads along and in the 
order of the direct connections in the interaction network (causality graph). In 
terms of example in Sec. 3, a failure of medical care would first affect the num- 
ber of infected persons, and in a second step the economic situation, public 
transport, and the disposal of waste. 

For a description of the catastrophe dynamics, let us assume that Pi{r) de- 
notes the impact on factor i at time r and Wji the rate at which this impact 
spreads to factor j, while is a damping rate describing the mitigation of 
the catastrophic impact on factor i by disaster management. In this case, it is 
reasonable to assume the dynamics 

dP - - 

— = (W - D)P(r) = LP(r) (30) 

with D = {S,jD,), L = (L,j) = {Wij-6ijDi), and P(r) = (Pi(r)). The symbol 
5ij represents the Kronecker function, i.e. it is 1 for i = j and otherwise 0. 
When no better information is available, we may assume that the spreading 
rate Wij is proportional to the strength |My | of the direct influence of factor 
j on factor i. With a constant proportionality factor c this means 

Wij^c\Mij\. (31) 

The formal solution of equation (30) for a time-independent matrix L is given 

by 

oo k 

P{t) = exp(Lr)P(0) = ^ v^L^P(O) = B(r)P(0) . (32) 

A;=0 

That is, B(r) describes the spreading of an event in the causality network 
(interaction network) in the course of time r, while P(0) reflects the initial 
impact of a catastrophic event. A random series of catastrophic events can be 
described by adding a random variable ^{t) to the right-hand side of Eq. (30). 

When we assume 

Di^^W,,, (33) 
j 

equation (30) is related to the Liouville representation of the discrete master 
equation. In this case, we can apply all the solution methods developed for 
it. This includes the so-called path integral solution [23], which allows one to 
calculate the occurence probability of specific spreading paths. This has some 
interesting implications. For example, the danger that the impact on sector io 
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affects tfie sectors ii,i2, ■ ■ ■ ,in in tfie indicated order is quantified by 

(34) 

Moreover, the average time at wfiidi tfiis series of events fias occured can be 
calculated as 

n ^ 

T{io ^ ii ^ ' ^n) = X] ' ^^^^ 

1=0 

and the variance of this time is determined by 

n ^ 

e(io ^ii^ >i^)^Y, ■ (36) 

z=o \^h) 

That is, Eq. (30) does not only allow to assess the likelyhood of certain series 
of events rather accurately, but also their approximate appearance times. In 
other words, we have a detailed picture of the potential catastrophic scenarios 
and of their time evolution, which allows for a specific preparation and disaster 
management. 

In the following, we do not want to restrict to the case (33). If 

A < E (37) 
3 

for all i, the damping is weak and the solutions Piir) are expected to grow 
more or less exponentially in the course of time, which describes a scenario 
where control is lost and the catastrophe spreads all over the system. In many 
cases, we will have 

Di>Y,Wji (38) 

for all i, i.e. the impact of the catastrophe on the system decays in the course of 
time, and \\m.-,-^Q Piij) 0. This determines, how strong the damping effects, 
i.e. the means counteracting the catastrophe have to be chosen. In terms of 
Sec. 7, this concerns the specification of the parameters Vn^. 

Finally, it may also happen that > J2j for some factors i, but < 
Y^j Wji for others. In such situations, everything depends on the initial impact 

P(0) and on the matrix B(r). However, in all these cases, Eqs. (34) to (36) 
remain valid. 



5 An Excitable Media Model of Disaster Management 



The damping effects Di are, to a large extent, related to the forces counter- 
acting the catastrophe. Therefore, we will now develop a dynamical model for 
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these, while our previous considerations assumed some more or less constant 
value of Di. Let us denote by the quantity of human forces (e.g. police, 
fire fighers, or military) ready for action, or the quantity of materials (e.g. 
technical or medical equipment) ready for use to fight the catastrophe. The 
index k distinguishes different kinds of forces or required materials. We will 
assume the following equation: 

^ = ^f- ± A,Ar,±e-^'=^ - \PmkNu{T)A,\T) . (39) 

The first term on the right-hand side describes the quantity of (human) 
forces and materials, which were exhausted or not usable, but become available 
again after an average time period of T^. The second term delineates reserve 
forces of quantity Nf, which are activated from the "standby mode" at a rate 
Afc after the occurence of a catastrophe (plus sign), while they are removed 
after the recovery from the disaster (minus sign). In most cases, A^^ < A^^, due 
to possible fatalities. The third term describes the activation of the forces k to 
fight the problems with factor or sector i. For simplicity, it is here assumed to 
be proportional to the strenth |Pj| of the catastrophic impact on factor i, with 
proportionality factors Vik which refiect the priorities in disaster management 
and the speed with which the forces or materials k become available for i. 
The exponent / allows to distinguish different cases: When the impact of the 
catastrophe on factor i is known, we may assume I — 0. However, when the 
active forces are assumed to order more forces, an exponent I > can make 
sense as well. 

The quantities of forces in action or materials in use change in time ac- 
cording to the equation 



(40) 



The first term on the right-hand side is due to the available forces /c, which 
are activated for fighting the catastrophe, while the second term describes 
a reduction of the active forces, as they become exhausted or damaged and 
require rest or repair after an average time period of T^. The third term 
describes unrecoverable losses such as fatalities or unrecoverable damage of 
materials, which are assumed to occur with a rate iPjli^ifc proportional to the 
catastrophic impact |Pj|. 

The quantity of exhausted or damaged forces is described by the following 
differential equation: 

dRk _ Ak{T) _ Rkjr) 

dr ~ T« ■ ^ > 

Herein, is the rate at which the forces k become exhausted or damaged, 

while is the recovery or repair rate. 



14 



The above model already contain many effects which are typically relevant in 
practical situations. We should note that it is related to models of excitable me- 
dia developed to describe chemical waves, the propagation of electrical pulses 
in heart tissue, LaOla waves in human crowds in stadia [20], or the spread- 
ing of forest fires [4]. These models typically contain three different states: 
an excitable one, an active one. and a refractory one. In our model of disas- 
ter management, refractory states are described by the variables Rk, active 
states by the variables A^, and excitable states by the variables N^. Due to 
this analogy, we expect to find certain pattern formation phenomena for our 
model of disaster management. In a forth-coming paper, this aspect shall be 
investigated in more detail. 



6 Some Analytical Results 



In the following, we will try to get an idea of the possible behavior of the 
excitable media model of disaster management suggested in Sec. 5. The sta- 
tionary state of this model is, for i/j^ = 0, given by 

^ = ^ = E \Pi\VikNkA,' = const . (42) 

In order to investigate the sensitivity with respect to small perturbations, we 
will carry out a linear stability analysis of the simplified model with one sector 
i, \Pi\ ~ const., = 0, and one kind k of forces. Dropping the subscripts and 
defining P — |Pj|, the resulting coupled set of differential equations is: 



f = ^-™. (43) 

'^ = PVNA--^, (44) 
dR _ A _ R 

^ - ^ - ^ ■ ^45) 

Its stationary solution is given by N{t) = Nq, R{t) = Rq, and A{t) = Aq 
with 

^ = ^ = PV'/V„V. (46) 

It is stable with respect to disturbances, if all eigenvalues A are non-positive. 
These eigenvalues can be calculated in the usual way. Assuming 

N{t) = No+5N e^^ , A{t) = Ao+5A e^^ , and R(t) = Ro+5R e^^ , (47) 
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we find the following eigenvalue problem for the amplitudes SN, 5A, and 5R 
of the deviations from the stationary values A^O) and Rq: 



A 



^6N^ 

5A 
y5Rj 



-PVAr^ 



V 



-PVNolAo^-' 1/T« 
PVAo^ PVNolAo^-^ 




SA 



The eigenvalues A are the solutions of the characteristic equation 



^I'-^-A 



i-PVAo' - A) I^PVNolAo'-' - ^ - 

+ P'V'NolAo'^-' (-J^-a)=0. 



PVAq^ 



The three solutions are 



(48) 



(49) 



A 



1/2 



PVAo'-\lNo-Ao) 



PVAo^-'{lNo-Ao) + ^-^ 



PVAo^ 



(50) 



and A3 = 0. Taking into account Eq. (46), which imphes PVAq'' ^ — 
l/{NoT^), this becomes 



A 



1/2 



2T^ V A^o 



2TR \ 



1 



2T^ V A^o 



2TR 



An 



(51) 

A detailed analysis of this expression shows the following: The system behaves 
unstable with respect to perturbations, when the real part of one of the above 
solutions becomes positive, which is the case for 



(52) 



Otherwise (apart from the case of marginal stability resulting for the equality 
sign), perturbations are damped, but one can distinguish two subcases: For 




</_l_^<niin|— ^, _, -_ + 2 



A^ 
No 




, (53) 



the resulting solution is complex, corresponding to damped oscillations, while 
the system behaves overdamped in the remaining case, where perturbations 
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fade away without any oscillations. For disaster management, the linearly un- 
stable case and the case of damped oscillations are both unfavourable. There- 
fore, / should be small enough. Otherwise, if active forces recruit other forces, 
the resulting "autocatalytic effect" may cause instabihties or overreactions 
in the supply with forces and materials. This effect is most likely for disas- 
ters which nobody was prepared for, where the recruiting mechanism plays 
the most signifant role. It may explain the suboptimal distribution of forces 
observed in these situations [24]. 



7 Connection with Supply Networks and Production Systems 



Finally, we have to specify the influence of disaster mangement activities on 
the damping D-i of the impact Pj, which a catastrophic event has on sector i. 
Let us assume that we have K different kinds of forces, materials or technical 
equipment. With k & K, we will indicate that the forces or materials k can 
substitute each other, i.e. K summarizes equivalent forces or materials. On 
the other hand, certain actions require the simultaneous presence of different 
supplementary kinds of forces and materials. Let us assume that the quantities 
simultaneously required to reduce the problems with factor i are represented 
by the coefficients CiK- The units of shall be chosen in a way that the 
following equation holds: 

AM = (l-iOmin(5«LMi5WMl)l (54) 

where VikNk{T)Ak{T) is the rate of activating forces k to mitigate the situa- 
tion of factor i. This equation reflects that, if only one of the required forces or 
materials is missing, no successful action can be taken. Moreover, there may be 
a loss Li of efficiency, e.g. due to queueing or limited capacities (0 < Lj < 1). 
The formula is analogous to that for production systems, where a product 
cannot be finished, as long as some required part or worker is missing, and 
where finite storage capacities may cause losses [21]. Therefore, this formula 
delineates the inefficiencies in disaster management, which occur when forces 
and materials are distributed in the wrong way. It has sometimes been re- 
ported, that too many forces have been located at some place, and missing at 
others [24]. Here, models designed to optimize supply networks could help to 
optimize the efficiency of disaster management [22]. 

Note that, for disaster management, a slight generalization of formula (54) is 
in place, as improvization may cope with a lack of certain materials or forces. 
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It is reasonable to assume a generalized function Gq with 

A(r) = (l-L.)G,({5-^^l|^Ml)Ml)j) (55) 

and 



Herein, the minimum reflects the worst case, while the sum over K describes 
the best case (if we neglect non-linearities, which may sometimes arise due to 
synergy effects). The specification 




T.k^K\P^\V^kNk{r)Ak\r)\\ \^ (i:keK\Pi\VikNk{r)Aj}{T 
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CiK J / V K \ J 

' (57) 

describes both extreme cases. The sum over K corresponds to g = 1, while 
the minimum results for q —oo. Hence, a variation of the parameter q al- 
lows to investigate different possible scenarios lying between the best case and 
the worst case. In order to model synergy effects between different forces, one 
would have to add nonlinear terms, e.g. bilinear ones. However, this would in- 
troduce a large number of additional parameters, which are even harder to es- 
timate than the first-order effects included in our model. Note that we already 
have non-linearities in our model, namely the products \Pi{T)\Nk{T)Ak{T) and 
the function Gg. 

The framework of supply networks also allows one to move from the semi- 
quantitative description of disaster management in Sees. 2 to 4 to a fully 
quantitative one, if the required data are available (while we work with as- 
sumption (31) otherwise): One simple case of the supply network model pro- 
posed in Refs. [21,22,25,26] corresponds to the dynamic input-output model 

-J^-I:{S^j-cM{r) (58) 

with 

dQj^ ^ Vj{Nj)-Qj{t) ^^g^ 
dr Tj 

where iVj denotes the inventory (stock level) of factor or product i, Vj{Nj) 
the desired and Qj > the actual throughput of sector j, Tj the adaptation 
time, and Cij > the quantity of factor i needed per throughput cycle. (For 
details see Ref. [22].) One could also say, Cij are the entries of the input-output 
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matrix mcasTircd in economics, and CijQj is the flow of tlie quantity generated 
by factor i to factor j. In the hmit of short adaptation times Tj ^ 0, the above 
equations reduce to 

-^-T.i5,,-c,,)V,{N,). (60) 
j 

Let us assume that the stationary state of this supply system is given by 
Ni{t) = Nf. Moreover let us denote the deviations from the stationary state 
by 5Ni{t) = Ni{t) - 7V°. With 

V,{N,) ^ V,{N^) + 6N, = A, - B, SN, , (61) 

the hnearized version of Eq. (60) reads 

dSNi 



dr 



J2{Wij-BjS,j)SNj{T), (62) 



where Aj = Vj{N^), Bj = -dVj{Nf)/dNj > 0, and Wij = CijBj. Here, we 
have used that, for the stationary solution N^, 



^{S,j-c,,){A,-B,N^)^0. (63) 
j 

In Eq. (62), we have Bi = J2j ^ji because of J2j Cji — 1- Taking into account 
the additional contribution (55), we finally obtain the set of linear equations 

dp 

5^ = (W - D)P(t) = LP(t) (64) 
ar 

with W = {Wi,), D = {5rj{B, + A)), L = W - D, and P{r) = {P,{t)) = 
(SNi^r)). Although Pj can become negative due to some catastrophic impact, 
Eqs. (35) and (36) for the average occurence times and their variance still 
remain meaningful, when Di is replaced by {Bi + Dj). Hence, they can be 
used to estimate the time of impact on other factors or sectors in the supply 
network. 



One particularly important aspect of supply networks is their sensitivity or 
robustness with respect to perturbations. It is, for example, known that sup- 
ply chains may suffer from the so-called bullwhip effect, i.e. small temporal 
variations in the demand may cause large variations in the supply. This insta- 
bility leads both to undesireable delays in delivery at some places and large 
stock levels at others [21,22,27]. In disaster management, this effect can have 
serious consequences. However, anticipation is known to efficiently stabilize 
the dynamics of supply networks [25,26]. As the formulas from Sec. 4 can be 
used to estimate the approximate time at which certain factors are likely to 
be affected, they can help to optimize the supply chain management, in par- 
ticular to stabilize the supply of forces and materials in time. It is certainly 
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reasonable to have forces available in time to fight the spreading of the dis- 
aster to other sectors, rather than sending them all to the places which are 
already devastated. The philosophy is to reach an anticipative disaster man- 
agement rather than having a responsive one. To model anticipation, the term 
T,iPi{T)VikNk{T)Ak\T) has to be replaced by 

\Pi{T + ATMkNk{T)Ak\T) , (65) 

i 

where At denotes the anticipation time horizon. 

Even more interesting is the robustness of supply networks with respect to 
structural changes, e.g. when some supplier fails to work or to deliver. This may 
be investigated by changing the coefficients Vik and CiK, which characterize 
the supply network. The influence of the topology of the supply network on 
its robustness, reliability, and dynamics is presently under investigation [26]. 
It has, for example, been noticed that supply ladders are more robust than 
linear supply chains or supply hierarchies (see Fig. 2). It is not surprising that 
the redundance of supply ladders, i.e. the availability of alternative delivery 
channels, stabilizes the system compared to a linear supply chain. We note, 
however, that hierarchical systems are very common in disaster management, 
and better alternatives are expected to be found in a research project that we 
presently pursue. 

Some insights regarding the robustness of networks have already been gained 
for the world wide web and other networks [14]. However, it is questionable 
whether the results for small-world, scale- free or random networks can directly 
be transfered to disaster management. Further research in this direction would 
be very helpful. 



8 Summciry and Outlook 

In the past, physics has made significant contributions to the understanding of 
catastrophes. This concerns, for example, the statistics of extreme events and 
avalanche effects. In this paper, we have tried to indicate how physics could 
contribute to disaster management [28]. We have sketched a rather general, 
semi-quantitative approach for the assessment of interaction networks which 
can serve for decision support, as it allows to compare potential optimization 
measures including their side effects. The method takes into account feedback 
loops and can be generalized to a dynamical model, which is suitable to esti- 
mate likely sequences of events and the times at which they are expected to 
materialize. We have also described the dynamics of disaster management in 
a way similar to excitable media, distinguishing an excitable state ( "ready to 
go"), an active state, and a refractory (exhausted) state requiring recovery. 
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Figure 2. Illustration of different supply networks: (a) linear supply chain, (b) "sup- 
ply ladder" , and (c) hierarchical supply tree. The supply ladder is particularly robust 
because of its redundant links and nodes. 

As successful actions need the simultaneous presence and/or action of several 
specialized forces and particular equipment, there was also a direct connec- 
tion with the management of supply networks. We have pointed out that the 
management of disasters and supply chains can be considerably improved by 
anticipation, based on the formulas in this paper. Moreover, the robustness 
[29] crucially depends on the structure of the supply network. In an on-going 
study, we investigate the optimal network structure to achieve robustness with 
respect to dynamical and structural perturbations [26]. The statistical physics 
of networks and graphs [14] is expected to make significant contributions to 
this. Some relevant aspects for the optimization of organizations and work 
groups have already been studied [30]. 

Our proposed approach connects to several methods and fields from statistical 
physics, such as the master equation [16,23], excitable media [20], and the dy- 
namics of transport (supply) processes [22] . It is also in the tradition of system 
dynamics [19], which has, with some success, been used to anticipate future 
problems of society [31]. In such kinds of studies, it is reasonable to carry out 
a sensitivity analysis [32] and to investigate the impact of random effects [16], 
which do, of course, play a significant role for the dynamics of catastrophes. 
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Apart from stochastic methods, one may in the future also apply elements of 
fuzzy logic [33] in order to describe the vague knowledge and soft facts, on 
which disaster management is often based. Insufficient, inconsistent, and un- 
certain information is one of the typical complications of disaster management, 
which makes it difficult to assess alternatives and to take the best decision, 
in particular under often very tense time constraints. In the future, informa- 
tion theory [34] is expected to make some valuable contributions to the design 
of decision support systems which can integrate inconsistent information and 
handle incomplete information [35]. 

In summary, statistical physics offers various promising concepts to develop 
and improve methods of disaster management. We think that the theory of 
self-organization [36] is particularly promising for this, having in mind prin- 
ciples such as synchronization [37], distributed control [38], and optimal self- 
organization [39] . It is expected that an application of these principles would 
lead to a more flexible, efficient, and robust disaster management compared 
to the present centralized or hierarchical concepts. 
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